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Abstract 

This paper, the last in the series, continues developing the nonlinear constitutive relations for non- 
isothermal, compressible, solid viscoelasticity. We initially discuss a single integral approach, more 
suitable for the glassy state of rubber-like materials, with basic functionals involved in the thermodynamic 
description for this type of viscoelasticity. Then we switch our attention to analyzing stability constraints, 
imposed on the general formulation of the nonlinear theory of solid viscoelasticity. Finally, we discuss 
specific (known from the literature or new) expressions for material functions that are involved in the 
constitutive formulations of both the rubber-like and glassy-like, complementary parts of the theory. 
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1. Introduction 

A thermodynamic theory for linear viscoelastic phenomena that accounts for thermo- 
rheological complexity has been developed in the first paper of this series (Freed and 
Leonov, 2002i). There, different temperature dependences for various parts of the 
relaxation spectra were modeled differently for the rubbery and glassy types of 
mechanical behavior for solid polymers. It has been demonstrated experimentally and 
discussed in various papers that matching of both time-temperature scales is quite 
possible. 

In the second paper of the series (Freed and Leonov 20022), a theory of solid 
viscoelasticity for large strain analysis has been developed. This nonlinear theory, 
originated from linear viscoelasticity with a discrete relaxation spectrum, is suitable for 
analyzing the non-isothermal and (generally) compressible behavior of rubber-like 
materials with large strains. 

For the sake of convenience, we now briefly present the main results from our 
previous paper (Freed and Leonov 20022) for nonlinear viscoelasticity with a single (bulk 
and shearing) nonlinear relaxation mode. The thermodynamic state variables in this case 

are: the temperature T , the density p , the “relaxing density”, p ' , the “incompressible” 

part of the total Finger tensor, B , and the incompressible part of the relaxing Finger 

A _ 

elastic strain, b . Completely bulk/shear decoupled, nonlinear, non-isothermal CE’s for 
solid rubbery viscoelasticity with a single relaxation mode, derived previously (Freed and 
Leonov, 20022) consist of: 

(i) A completely decoupled free-energy function, 


F = F? (T, p) + f; (T, p‘ ) + F* (T, /',/*) + /£ (T, /, , / 2 ) ; 


(ii) A corresponding, additive presentation of the stress, 


£ = -(Po + Pi)£ + £ 0 + £ 1 ’ Po = p 2 dF I dp, p x = pp‘dFldp‘ , 
£ 0 = 2p[F?(£-l/3lf£>-F?(B'-l/3I 2 B S)], 
a x = 2 p[F? (b - 1 / 37 , £) - F± (b~' - 1 / 31 2 £)] ; 


(iii) Continuity and (shear) kinematical relations, 


( 1 ) 


( 2 ) 


V 0 

dkip/dt = -tre , B = B-B-e-e-B = 0(e = e-l/3Stre); (3) 


(iv) Bulk and shear evolution (relaxation) equations, 


d\np i /dt + [«(/>')/ &l (T)]pp f dF^ / dp\T, p‘) = d\np/dt. 


(4a) 



(4b) 


b-b-e-e-b+ \l2[p{i x J 2 )ie x {T)]-[^ -1/3(7, -I 2 )b-S] = 0. 

These equations have been shown (Freed and Leonov, 20022) to easily extended to the 
multimode case with an arbitrary number of nonlinear relaxation modes. Additionally, 
there is a strain-induced anisotropy in the thermal conductivity, which in the simplified 
version depends on temperature T and the actual component B of Finger strain. 

In the rubbery region, the well-known, WLF time-temperature superposition 
principle, originated from entropy elasticity and uniform temperature scaling for 
relaxation modes, is applicable for scaling of the linear viscoelastic spectrum. It is also 
applicable in the nonlinear region of deformations, too. However, in the glassy state there 
is another temperature scaling related to the possible fractal properties of polymers in the 
glassy region. Therefore, viscoelastic effects in the glassy region are well described by 
die stretched exponential and/or the fractional derivative approaches. Using a discrete 
relaxation spectrum approach in this region is awkward, as it usually includes a lot of 
relaxation modes. In order to expand this type of formalism into the nonlinear region of 
deformation, a slightly modified single integral approach to constitutive modeling of the 
K-BKZ type (Kaye, 1962; Berstain, Kersley, Zapas, 1963) is introduced in the next 
Section of this paper. This approach has a long history of use for describing nonlinear 
viscoelastic properties of solid polymers. Some new approaches have also been recently 
attempted. They combine continuum mechanics with more physical insight in the 
behavior of materials in the glassy region. For instance, Drozdov (2000i) recently 
proposed CE’s for nonlinear viscoelasticity of glassy polymers based on the concept of 
traps. He also developed a model of cooperative relaxation in finite viscoelasticity for 
amorphous polymers (Drozdov, 2OOO2). Huang (1999) recently extended Drozdov’s 
approach to the non-isothermal case. Another approach that has also been developed 
recently, is based on non-equilibrium thermodynamics and the concept of the soliton-like 
propagations of the specific structures in glassy polymers (Cho and Kim, 2000). 
However, these approaches contain a lot of fitting parameters, many assumptions, and 
have been tested only for the one deformation (ID) case, as uniaxial extension. 

Another issue addressed in this paper deals with stability constraints, which 
secure fundamental properties of CE’s, such as obeying the second law of 
thermodynamics, and well-posedness. We will also discuss some recent results in this 
area in the third Section of this paper. 

At the end of this paper we discuss formulae proposed by various authors for the 
purpose of specifying material functions and parameters that appear in our viscoelastic 
theory. 

2. Single integral constitutive equations for viscoelastic solids 

The general CE for viscoelastic materials, presented in a form of hereditary functionals 
with fading memory, mostly applicable to viscoelastic liquids, has been established in the 
field of rational mechanics and thermodynamics (Truesdell and Noll, 1992). In order to 
apply to a viscoelastic solids the fading memory concept using a single integral CE, it is 
reasonable represent the Helmholtz’ free energy as a sum of the equilibrium free energy 
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depending on temperature and actual strain, F°(T,B ) , and the liquid-like free energy in 
a form of a single integral hereditary functional with the integrand depending on 
temperature and a relative total strain, F l {T,BJ . We assume further that both 

components of the free energy vanish in the unloaded state, where B = B r =6. For 
isotropic materials, the free energy can then be expressed as: 

F = F\T-,If,i;,I s ‘)+ (5) 


Here the integrand in (5) depends on the “relative” Finger strain B^(t,z) which acts in 

the past time region { v } and degenerates into the unity tensor at the actual time r = / . 
Consequently, the kinematical equation for the evolution of the relative Finger strain 

B* holds here too, but with a specific initial condition: 

Ir '£“£ 'Ir = 2 ; IrU = £’ (6) 

The initial condition in (6) means that the strain, initiated at the current time t, 
“propagates” into the past. The stress tensor is defined now by the Finger formula but 
with a functional (Frechet) derivative: 


a = - = 2pB-dF°/dB+2p\B T -dF'/dBdT. (7) 

— - —oo 

Substituting into (7) the isotropic expressions for F° and F ] yields a BKZ type of 
formula, with an additional out-of-integral, equilibrium term: 


g = 2p[F l °B + F^I l B B-B 2 ) + F 3 °I 3 B S]+2p jdr [i* 1 * +F}(I*B r -g t )+F}l 3 B T £] 


( 8 ) 


Fj = dF° /dl°; F; = F]{T-Il,t~r) = 8F l !dl* x ; pit) = p 0 (I 3 B y m ; p = p T (I 3 * T ) 


B \-l/2 


Here p T is the density at a past time r . To find the expression for the dissipation D we 
employ the expression for the entropy production under isothermal condition: 

D = TP s \ T =tr(a-e)-pdF/dt\ r (>0) 

Substituting here relation (5), taking into account the kinematical relation of (6) yields: 
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(9) 


D = -p jA(T-J? t X,I 3 B T ;t-T)dT 


Here Ff means a differentiation of the function Z 7 ' taken with respect to the third (“pure 
time”) argument. Note that the positive definiteness of the multipliers in the integrand of 
(9) necessitates the inequality /& < 0 . This is the nonlinear extension of the decaying 
exponential factors in the Prony series, well known in linear viscoelasticity. 

Instead of invariants If and/.*, we now introduce density variables into the 
equilibrium and non-equilibrium components of the free energy as was done in the 
previous paper (Freed and Leonov, 20022): 


F = F\T,p! jP'(T;p/ 


( 10 ) 


Here the invariants Ij and/", are defined as in the incompressible case. Then the 
constitutive equation in (8) is presented as: 


t 

a = -/>£+<£ (tr<x = 0); p- p 1 dF°/dp + jdzdF ] /dp 

-ao 

a = 2p[F? ( J - 1 / 3/* £) - F° ( J"' -1/3 I*8j\ 

+2 p )dr [F?(B t -miZg-FXB' -1/3/' 0 ] . 


( 11 ) 


^ A ^ 1 

Here i 7 - and Fj are defined similarly to that shown in the previous paper (Freed and 
Leonov, 20022). The relations in (11) can be simplified, if there is an additivity (or 
uncoupling) assumption between the bulk and shearing responses of elasticity and 
viscoelasticity: 


F° =F?(T;p / + F l = Fj(T;pt p T ;t -t)+ F£(T;i t B T Jf T ;t -t) . (12) 


A. A. 

Note that the time dependencies of functions F 4 and/ 7 ), in (12) for the compressible 
(scalar) and “incompressible” (tensor) terms are in general different. 

The constitutive equation in (11) with assumption (12) can be further simplified 
when it is possible to employ a factorization called time-strain separability: 

Fl =m b (t-T)m P /Pr)> Fl=m^t-T)mi?J!r) (13) 

Here m h (t - r) and m sh (t - r) are respectively related to the bulk, K(t ) , and shearing, 
G(t ) , relaxation moduli defined in the linear theory of viscoelasticity as follows: 
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m b (t-z) = dK(t-z)l dz (>0), m sh (t-z) = dG(t-z)/ dz (>0) . 


(14) 


Thus the most simplified version of a single integral CE for viscoelastic solids has the 
form: 







dr* 



dr -p \dm*«-T)ldt f(T-J‘Jl)dr 


(> 0 ) 


cr = - pS^ + & (tr<r = 0) ; 


P = P 


2 


dF b / dp + ^m h (t- z)d(f> / dpdz 


(15) 

(16) 


g = 2p[/|”(|-l/3/,'£)-^(|" 1 -l/3/'£)] (17) 

+2 p r) Did, -1/3 i‘£)-/ 2 (£' -1/3 if,£)]dr (fj=df/d/? r ). 

—oo 


Due to the definitions given in (14), the dissipation in (16) is positive, because 

dm h (t-z)ldt = d 2 K(t -z) / dtdr (<0), dm xh (t -z)/8t = d 2 G(t - r) / dtdz (<0) . 

Comparing the single integral (SI) approach considered here with the differential 
approach analyzed in the previous paper (Freed and Leonov, 20022), we can conclude: 

(i) The SI CE’s are more rigid when compared to the CE’s of the differential type. As 
seen from equations (15) and (16), the dissipative function D is almost proportional to 
the integral part of the free energy. For the CE’s of differential type, these two 
fundamental functions are completely separated. This happens in part because of the use 
of relative strain measures in SI CE, instead of the hidden tensor variables in the 
differential CE’s. 

(ii) The viscoelastic SI CE’s are considerably “less nonlinear” than CE’s of the 
differential type. The latter can be, in principle, represented by a multiple integral 
approach. More nonlinearity, as e.g. a “strain clock”, could be introduced into the non- 
time separable versions of SI CE’s, however, at an expense of tractability. 

(iii) For viscoelastic liquids, SI CE’s with time-strain separability were proved to be 
unstable (Kwon and Leonov, 1995; Kwon and Co, 2001); however, for viscoelastyic 
solids, stable SI CE’s with time-strain separability can be constructed. Therefore, unlike 
the CE’s of differential type, it is impossible to use the SI CE’s in the continuous 
modeling of rubber viscoelastic properties, from liquid-like to solid-like via cure reaction, 
using the chemo-rheological approach (Mitra, 2000). 
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(iv) There are also some advantages of the SI CE’s: 1) relative simplicity since it is 
necessary to specify only one integrand (and out-of-integral, equilibrium term in the solid 
case); 2) the possibility to involve in the formulation a continuous viscoelastic spectrum 
with such features as fractional derivatives and “stretched exponents” (e.g. see Kolrausch, 
1847, 1854), valuable for glasses; and 3) the use of a rigorous approach to viscoelastic 
kinematics. 

Finally, it seems that a combination of the differential approach with a few modes 
and the SI CE’s for fast deformations may create a good opportunity for a unified 
description of nonlinear thermo-rheologically complex viscoelastic solids and liquids. In 
this unified scheme, the differential approach can take care of the rubbery region, with 
relatively slow deformations and potentially high strains, and with the WLF temperature 
scaling. The SI constitutive approach, with relatively small strains and another 
temperature scaling, is capable of describing the glassy region. 

3. Stability constraints for constitutive equations for viscoelastic solids 

Formulation of field equations for viscoelastic solids should obey some stability 
constraints. These constraints reflect the most fundamental properties of CE’s, related to 
the second law of thermodynamics, propagating various fast waves, and continuity along 
the time axis for any solution for a viscoelastic solid problem. We will show below that 
the stability constraints in the case of viscoelastic solids are exactly the same as in the 
pure elastic case, i.e. they are related only to Hadamard stability (or well posedness, or 
strong ellipticity, or strong hyperbolicity in case of wave propagation). 

In the theory of finite elasticity, the necessary and sufficient condition for strong 
ellipticity is the positive definiteness of the bi-quadratic fourth order form, 

A = A.. mn {E)x i y j x m y n > 0, ( 1 8) 


where the tensor A, Jmn is related to the elastic properties of the continuum. 

The rank of four tensor A,y m „ is symmetric in its first and second pairs of indices, 
and also by a transposition between these. For isotropic solids it has the form: 


dtju 

A iimn (B) = 2- 

ljmn ~ dH„„ 


8H„ 


dF 


dH, 


= 4 B, 


mk 


9 J 


dB 


kn 


p b « 


dF 

dB Uj 


(19) 


Here H is the Hencky strain measure. The condition (18) of strong ellipticity means the 
positive definiteness of the tensor A,j mn considered in terms of a bi-quadratic form. In the 
isotropic case, when F = F{T,I l J 2 ,I i ) , the stability constraint (18) was reduced in both 
incompressible (Zee and Sternberg, 1983) and compressible (Knowles and Sternberg, 
1977; Simpson and S.J. Spector, 1983) cases to corresponding sets of algebraic 
inequalities imposed on the first and second free energy derivatives with respect to the 
three basic invariants I k . In the case of the linear stability analysis of disturbances 
imposed on an arbitrary, nonlinear, and generally time dependent stress-strain 
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deformation field, condition (18) of strong ellipticity is changed to the condition of strong 
hyperbolicity: 


P&\Vi = A Jm n(iX k j Y Jn > 0. 


( 20 ) 


Here Q. = co-k- v is the frequency oscillations in the disturbance wave with aDoppler’s 
shift of the basic velocity field v = v(x,/) taken into account, k\ and v- are respectively 
the components of the wave vector and the vector of amplitude in the velocity 
disturbances. In the limit co -» °o, k -*• oo , one can consider disturbances as presented by 
monochromatic plane waves, being restricted to a local analysis of stability. It is also 
worth mentioning that the condition of stability (20) is exactly the same as that shown in 
(18). 

When a system is described by the quasi-linear set of PDE, the analysis of small 
disturbances has been shown to coincide with the more general method of characteristics 
(e.g., see Joseph, 1990). Positive definiteness of the forth order form in (20) means the 
possibility to continue at any time instant, any solution of the equations for finite 
elasticity in the positive direction of time axis. If this is impossible, the characteristics are 
turned back, and near a turning point, their derivatives go to infinity. Thus the Hadamard 
instability is a typical ill posedness accompanied by a blow-up instability. In the general 
case of compressibility, all possible types of waves of disturbances are allowed, including 
the longitudinal waves. However, when an elastic material is treated as incompressible, 
i.e. p = const , then the longitudinal waves of disturbances are forbidden because in this 
case the wave and the velocity amplitude vectors are orthogonal: £ y v y = 0 . 

When analyzing static stability, with infinitesimal (symmetric) strain disturbances 
£ imposed on a given finite elastic strain B Q as B = B q +s, one can obtain instead of 
the strong ellipticity condition (18), a condition of thermodynamic stability: 

A 7 - = Ajjm„ (B) £ jj£mn > 0. (21) 

Although the tensor Ay mn is still defined by (19), the condition of stability is different 
from that which is shown in (1 8). It is easy to see that the thermodynamic condition of 
stability (21) is only a necessary condition for the full condition of strong ellipticity (1 8). 
This is because the condition of strong ellipticity degenerates into the condition of 
thermodynamic stability for a particular case when the diadic matrix x,.y y in (18) is 
symmetric. In this case, instead of positive definiteness of the bi-quadratic form (18) for 
3D vectors, the thermodynamic stability condition reduces to the positive definiteness of 
the quadratic form (21) for 3D matrices. 

In the case of CE’s for viscoelastic liquids with instantaneous elasticity, the 
general stability constraints have also been revealed and analyzed (Kwon and Leonov, 
1995; Kwon, 1996; Kwon and Co, 2001). It was found that the Hadamard stability 
constraints for CE’s of viscoelastic liquids coincide with those found for the elastic 
solids. However, another source of instability, caused by dissipative terms in liquid-like 
viscoelastic CE’s, can generate dissipative instability even in the Hadamard stable case. 
Therefore the combined constraints have been imposed on liquid-like viscoelastic CE’s. 
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For CE’s of differential type the degree of freedom was high enough to allow a 
completely stable class of viscoelastic CE’s to exist. This class described well the various 
rheological properties for five polymers {Simhambhatla and Leonov, 1995; Leonov, 
1999) while obeying all the stability constraints. However, the combined stability 
constraints imposed on the single integral CE’s with strain-time separability are so stiff 
that a stable version of these does not exist (Kwon and Leonov, 1995). 

It is a great relief that the dissipative instability does not exist in solid-like 
viscoelastic CE’s. Therefore, only the Hadamard stability conditions of (18) or (20) must 
be satisfied to guarantee a stable constitutive behavior for viscoelastic solids, whenever 
the dissipative terms there subordinate to the second law. It should be mentioned that the 
same necessary and sufficient stability criterion of (18) or (20) for elastic terms in each 
mode of the multi-modal differential CE’s or in a single integral CE (with the strain-time 
separability admissible here) can be proven for the CE’s of viscoelastic solids. Strictly 
speaking, this means that the above Hadamard stability criterion is only sufficient. 
However, since the various viscoelastic modes and the equilibrium elastic term are 
almost independent, the sufficient stability criteria of (18) or (20) should be considered as 
close to the necessary one. 

The strong ellipticity (Hadamard stability) conditions are too awkward to apply 
them for newly proposed CE’s. Therefore several attempts have been undertaken to find 
non-trivial and simple enough sufficient conditions for strong ellipticity. We discuss 
below only two literature approaches to this problem. 

(i) Renardy (1985) derived the following sufficient condition of strong ellipticity for the 
incompressible viscoelastic Kaye-BKZ equation:. 

The CE’s for incompressible elastic and viscoelastic solids are strong elliptic, if the 
strain energy function F = F(T,I U I 2 ) , depending on two basic invariants, // and /?, 

can be presented as a monotonically increasing and convex function of ^/7j~ and . 
The condition for this convexity, 


d^F d 2 F f d 2 F T 

dfidtf ‘ dJT 2 dJT 2 ~ { dyftdypl ) ’ 


( 22 ) 


is reduced to the following: 

(F\\Fji- F X 2 )IiI 2 +\l2(F n F 2 I x + F 22 F X I 2 )+H 4F X F 2 > 0 
( F k ^dF/dI k ;F ik ^d 2 F/dl ,dl k ) 

Evidently, inequality (23) is less strong than the convexity condition, F n F 22 >F* 2 , 
imposed on F in terms of the basic strain invariants, I V I 2 . 

(ii) Kwon (1996) derived a sufficient condition of strong ellipticity for compressible 
viscoelastic liquids, which can also be used for both elastic and viscoelastic solids. Since 
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these conditions are considerably more sophisticated than Renardy’s , we refer the reader 
to the original Kwon paper (1996). 


4. Specification of material functions and parameters 


This Section discusses many formulae derived, or proposed at hoc, by various authors for 
material functions, which are involved in the general theory. 

4.1. Equilibrium (elastic) properties 


These properties in our simplified approach are described by two independent parts of the 
free energy function, the bulk part, F h depending on temperature T and density p , and 
the “shearing” part, F'* depending onT and two “incompressible” invariants, /, = trB 

and I 2 = trB ' (det B = 1 ). 

(i). Spenser and. Bauer (1946), and Spenser and Gilmore (1960) have established 
that the bulk properties of cross-linked rubbers can be well described by a van der Waals- 
like equation: 


(7» = ^ln| 

M 


p 0 (M/b-p) 


M 2 


(24) 


Here p is the density, T is the Kelvin temperature, M is the molecular weight, R is the 
universal gas constant, p 0 is the density at a reference state under temperature To. The 
parameters a and b in (24), treated as empirical constants, describe the effects of internal 
pressure and the self-volume of molecules, respectively. In accordance with relation (24), 
the dependence of pressure ponp and T is: 

p = p z dF b / dp = pRT /(M -b p) - a p 2 / M 2 (0 < p < M/b) . (25) 

The global (Hadamard) stability condition, F p + pF h pp / 2 > 0 , well known from 
thermodynamics of gases, reduces to the inequalities: 

p< p c =M /(3b), T>T c =(Sa)/(27bR), (26) 

Here p c and T c are the critical density and the temperature at the point of phase 
separation for gases. 

(ii) Consider now the specifications of potential F sh = F sh (T,I X ,I 2 ) for shearing 
equilibrium properties. To simplify the notations, we consider the incompressible case 
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below, where F sh is changed to F, and the modified invariants/, and / 2 are changed 

to /, and I 2 , i.e. F = F(T, /, , / 2 ) . Additionally we will use the notation: pF = W . 

The most popular mode is the Mooney-Rivlin potential: 

W = G/ 2[(1 - «)(/, - 3) + a(I 2 - 3)] , (27) 

where G and a are respectively the Hookean modulus and a numerical constant. The 
stability constraint provided by the incompressible stability theory (Kwon and Leonov, 

1995) is: 0 < a < 1 . However, when the Mooney potential is combined with an additional 
compressible part, the constraint on the parameter a is tougher: 1 / 3 < a < 1 (Kwon, 

1996) . The descriptive ability of this potential in the incompressible case has been widely 
discussed (e.g., see Treloar, 1975). 

To the authors’ knowledge, the widest range in descriptive capability of various 
experimental data for cross-linked rubbers has been achieved using the BST potential 
(Blats, Sharda, Tschoegl, 1974): 

W = (G 0 / 2 ri)I„ + GI m „ , /„ = nMl” - S ) . (28) 

Here Go and G are effective moduli, and n and m are numerical constants. It should be 
noted that using the generalized strain measure B" introduces an additional 
computational expensse, because the stress-strain relation with potential (28) can be 
effectively treated only in the principal axes of tensor B . This creates a computational 
burden when using FEM to solve a complex problem. The stability constraints for the 
potential (28) are unknown, mostly because of the generalized strain measure involved in 
its formulation. 

Simhambhatla and Leonov (1995) recently proposed another expression for 
elastic potential involved in their constitutive theory for viscoelastic liquids: 

w = ^[(1 - «)[(/, /3)" -1] + <*[(/, 13] r - 1]]. (29) 

The corresponding stress-strain relation is: 

a = -pS + G[(l - «)(/, / 3)"-' B-cc(I 2 / 3)"-' B~ l ] . (30) 

Here G is the Hookean modulus, while a and n are numerical constants. Expression (29) 
is a combination of the formulae used for Mooney-Rivlin and Knowles’ (1979) 
potentials, widely used for calculations of large deformations of elastic solids. Using the 
combined Hadamard and dissipative stability criteria, it has been established that the 
sufficient conditions for Hadamard stability of viscoelastic liquids using potential (29) 
are: n > 1, 0 < a < 1 . However, using the Renardy sufficient condition of strong 
ellipticity, results in softer stability criteria: n> 1 / 2 , 0 < a < 1 , suitable for elastic 

solids. 
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Several expressions for an elastic potential were recently proposed for very large 
deformations of elastic solids (e.g., see Arruda and Boyce, 1993; Bergstrom and Boyce, 
1998; Boyce and Arruda, 2000). These formulae were derived using an extension of the 
finite extensibility approach for polymeric chains (e.g. see Treloar, 1975). These 
formulae employ the inverse Langevin function, which is too complicated for solving 
complex problems. An easy semi-empirical way to involve these effects into 
consideration was recently proposed by Gent (1996). Gent was unaware that Warner 
(1972) had earlier proposed the same approach for approximating the inverse Langevin 
function. The Warner approximation is widely used for viscoelastic liquids as the 
“Warner spring” (e.g., see Larson, 1988). This type of potential and its corresponding 
stress-strain relation are: 


W = -(G / 2)(7. - 3) In ■ 
V 


i-/,//. 

1-3/7. 


cr = -p5 + G — 3/It -B. 
= = 1-7,//.= 


(31) 

(32) 


Here 7. (~100) is a limiting value for the first invariant 7i related to the finite extensibility 
of polymer chains. The potential in (31) is stable in the entire domain of admissible 
strains: 3 < 7j < 7. . It should be noted that when 7i approaches 7. , the concept of entropic 
elasticity is becoming invalid. Here die “solid-like” elasticity, related to the distortion of 
valence angles between monomeric units, needs to be brought into consideration. 

The potential of (31) might describe large elastic deformations very well. These 
usually occur in uniaxial extension or in situations close to that, which can be met in 
complex deformations with high-level concentrations of stresses. However, at relatively 
low strains, predictions of potential (31) contradict experimental data (Treloar, 1975). 
These data clearly indicate that elastic stress also depends on the second invariant / 2 . 
This means that in low and intermediate regions of elastic strain, the potential (29) might 
be more preferable. To combine both potentials, (29) and (31), we suggest an 
interpolation formula for the elastic potential and a corresponding expression for the 
stress-strain relation of the type: 



- (1 - a)[(7. / 3)" - l]ln 

//.)■] 

+«[(/,/ 3)” - 1] 

In 


,i-(3//.r J 



^_^ + 3 M .[(7./3)"-l]- ( 7,/7.r 
= = 7. !-(/,//.)" = 2 


(33) 

(34) . 


Here G is the Hookean modulus, and 7* (-100), n and a are three numerical parameters. 
It is seen that the formulae (33) and (34) practically coincide with (29) and (30) at 
relatively small strains, but at large value of 7, , only the first, singular terms, in (33) and 
(34) are dominant. When n = 1, formulae (33) and (34) reduce to a combination of 
Warren-Gent and Mooney-Rivlin potential. Using the Renardy stability constraints, it is 
easy to establish that in the whole admissible wedge-like domain of invariants. 
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3 < /, < I,, 1 2 > 3 , when I, » 3 , these sufficient stability constraints are almost the 
same as for the potential in (29), i.e. 0 < a < 1 , n > 1/2. 

42. Non-equilibrium (viscoelastic) properties 

These properties include the linear viscoelastic spectrum, and the stress-elastic strain 
( potential ) relations for both the differential and single integral parts of the combined CE. 

(i) The linear viscoelastic spectra, for bulk and shearing infinitesimal deformations can 
be established from a series of isothermal experiments (either dynamic, or creep, or 
relaxation) by using the Pade-Laplace technique (e.g. Simhambhatla and Leonov, 1993), 
after extracting the equilibrium (static) term. This method of treatment of experimental 
data can also reveal die thermo-rheological complexity, when present. Although some 
molecular models also exist, they usually do not describe well the relaxation spectra. 

(ii) For the large deformation and slow relaxing differential modes, both the nonlinear 
bulk and shearing viscoelastic properties are described by a quasi-equilibrium elastic 
potential and non-equilibrium relaxation functions in the respective evolution equations. 
A great deal of simplification has been achieved for incompressible CE’s for viscoelastic 
liquids, using “uniform modeling” (Simhambhatla and Leonov, 1995; Leonov, 1999). 
Uniformity means that the same specifications apply to all nonlinear relaxation modes, 
with the same nonlinear numerical parameters, while keeping the linear viscoelastic 
parameters established in the linear viscoelastic region. In addition, uniformity provides 
nonlinear viscoelastic CE’s with easily formulated stability constraints. This uniform and 
stable modeling has been successfully tested with various modes of deformation for five 
common polymeric melts (Simhambhatla and Leonov, 1995; Leonov, 1999). Moreover, 
in this type of modeling, a scaling approach introduced by Leonov and Padovan (1999) 
can also be used, that provides a great computational economy when solving problems 
with complex deformations. 

In the incompressible case (or the shearing part of CE’s), two material functions 
have to be specified in the A-th nonlinear relaxation mode: the mode’s firee-energy 

function F k , scaled by the Hookean elastic modulus G k ~ pRTg k , and a determining 

scalar J3 k ; both functions depend on the basic invariants l[ k) and I 2 k) . A good descriptive 

capability is achieved here when F k is used in the form of (29), where G k is used instead 

of G , and l[ k) and I 2 k) instead of /, and I 2 . Uniformity here means that the values of 
the nonlinear numerical parameters a and n are mode independent. Uniformity in 
modeling functions J3 k means that J3 k - f(I k ,I k ), with the same, mode independent 
numerical parameters characterizing this function. The existence of a linear viscoelastic 
limit imposes the additional condition: J3 k -» 1 as I } -> 3 and I 2 -» 3 . In many cases, 

the simplest “power law” specification, J3 k = (/* / 3) m , where m is a mode independent 
numerical parameter, serves very well for modeling. 
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(iii) As mentioned, relatively low strain and fast-relaxing processes could be formally 
described by the CE of single-integral type (15)-(17). Here, different approaches to the 
linear viscoelastic response can be used, including fractional derivatives and fractal 
arguments (see the book edited by Hilfer, 2000). It is also possible to use here a more 
compact description of the linear relaxation spectrum, such as the Kolrausch (1847,1854) 
relaxation function: 

G(t) = G 0 exp[-(/ / 1 0 )'] , (3 5) 

with a positive numerical parameter s . 

4.3. Strain induced anisotropy in thermal conductivity 

As shown by Freed and Leonov (20022), our thermodynamic arguments resulted in a 
prediction of strain-induced anisotropy in the heat conductivity at large total 
deformations, i.e. in the occurrence of isotropic tensor dependence rc(B) . Some 

additional physical arguments (Hellwege, Henning and Knappe, 1963; Novichionok and 
Shulman, 1971) resulted in the simple expression for components of heat conductivity in 
the uniaxial extension: 

l/k ? +2/k ± =const = 3/tc 0 . (36) 

Here k 0 is the heat conductivity coefficient in isotropic (unloaded) state, kp and k j. being 

respectively the heat conductivities in the longitudinal and cross-sectional directions. 
Formula (36) has been experimentally confirmed for various elastomers. 

To describe (36) with the use of a general isotropic tensor function k = fc(B) we 
neglect compressibility effects at large deformations, and assume a quasi-linear 
dependence of the form: 

k = tc 0 <pB . (37) 

Here the unknown scalar factor tp can depend on temperature and the basic invariants /, 
and I 2 of tensor B . In uniaxial extension, the left-hand side in (36) is represented as 

trig' . Then the relations (36) and (37) yield: 

l/*r p + 2/*r x = trig' =^L = A_ = A. 

K 0 <p K 0 <p k 0 

This means that (p = I 2 l 3 , and the final expression for the heat conductivity tensor is: 

g(B) = l/3K 0 I 2 B. (38) 
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This remarkable simple result needs to be tested experimentally for other deformations. 


6. Conclusion to the series 

1. Based on the quasi-linear approach of non-equilibrium thermodynamics, a general 
nonlinear theory has been developed for arbitrary non-isothermal deformations of a 
viscoelastic solid with one relaxation mechanism. 

2. The theory has then been simplified using several modeling arguments: 

(i) Decoupling of the free energy into a sum of equilibrium and quasi-equilibrium 
parts; 

(ii) Decoupling each component into a sum of compressible and incompressible 
parts; 

(iii) Decoupling of the evolution equation for the hidden tensor variable (elastic 
strain in the viscoelastic mode) into two, almost independent evolution 
(relaxation) equations for bulk and “shear” components; 

(iv) Modeling of anisotropic heat conductivity as induced by the equilibrium 
(total) strain. 

It should be mentioned that the decoupling of bulk and shear phenomena is usually 
treated as evidence that the bulk (compressibility) effects are small. This is definitely true 
for large elastic deformations of rubber-like materials, but incorrect for small 
deformations, when the shear and bulk effects happen to be of the same order. 
Nevertheless, at small deformations, where the linear theory is approximately valid, 
additivity of the bulk and shear deformations remains correct. Therefore the proposed 
decoupling of bulk and shear effects seems to be approximately valid over the total range 
of deformations. 

3. The models in 2 have been extended for a nonlinear decoupled multi-mode CE’s. 

4. Some thermodynamic arguments have been used for a solid-like Kaye-BKZ type of 
single integral equation with out-of-integral equilibrium term, when decoupling 
equilibrium and non-equilibrium properties. The time-strain separable version of the 
single integral approach has also been discussed. 

5. It was suggested to combine the model in 3 with the single-integral approach for the 
description of fast, weakly nonlinear relaxation processes in the glassy polymer region. 

6. Stability constraints have also been discussed for the proposed CE’s. It has been shown 
that unlike the CE’s for viscoelastic liquids where the strong ellipticity (Hadamard ) 
stability constraints should be complemented by some criteria for dissipative stability, 
only strong ellipticity criteria should be used for stability of CE’s for viscoelastic solids. 
Special attention was attached to the sufficient stability criterion. 

7. Several formulae, some of them being new, have been proposed and discussed to 
specify the equilibrium and non-equilibrium terms in the solid viscoelastic CE’s, 
including the derivation of a simple formula for strain induced anisotropy in heat 
conductivity. 

Finally, predictions of this theory for nonlinear solid viscoelaticity, developed in 
this series of papers, need to be tested experimentally for several well defined and easily 
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checked problems. The work on this is in progress, and results will be published 

elsewhere. 
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